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In this paper, we consider the problem of computing the self-force and self-energy for a static scalar 
charge in a wormhole space-time with throat profile r(p) — \l p 2 + a 2 for arbitrary coupling of the 
field to the curvature. This calculation has previously been considered numerically by Bezerra and 
Khusnutdinov [T], while analytic results have been obtained in the special cases of minimal (£ = 0) 
coupling [2] and conformal coupling 1 (£ = 1/8 in three dimensions). We present here a closed form 
expression for the static Green's function for arbitrary coupling and hence we obtain an analytic 
expression for the self-force. The self-force depends crucially on the coupling of the field to the 
curvature of the space-time and hence it is useful to determine the dependence explicitly. The 
numerical computation can identify some qualitative aspects of this dependence such as the change 
in the sign of the force as it passes through the conformally coupled value, as well as the fact that 
the self-force diverges for £ = 1/2. From the closed form expression, it is straight-forward to see 
that there is an infinite set of values of the coupling constant for which the self-force diverges, but 
we also see that there is an infinite set of values for which the self-force vanishes. 



I. INTRODUCTION 

Wormholes are topological bridges connecting differ- 
ent universes or distant regions of the same universe. In- 
terest in wormhole space-times dates back to 1916 [3J, 
prc-dating interest in black hole space-times. Its modern 
popularity is owed primarily to the work of Morris and 
Thorne [3] who investigated the idea of using so-called 
"traversable wormholes" as a means for time travel. Mor- 
ris and Thorne showed that such space-times require a 
stress-energy tensor that violates the null energy condi- 
tion, that is they require the existence of exotic mat- 
ter. Wormholes have subsequently played a central role 
in the investigation of causality violation and the status 
of energy conditions as physical laws of nature. Com- 
prehensive reviews of wormhole physics can be found in 
Refs.[3E]. 

We will consider a particularly simple, ultra-static, 
spherically symmetric wormhole geometry described by 
the metric 



ds 2 = -dt 2 + dp 2 + r 2 (p)dn 2 , 



(1.1) 



where dQ 2 — d8 2 + sin 2 9 d<f> 2 is the line-element on the 
two-sphere S 2 and r(p) = \/ p 1 + a 2 is the profile of 
the wormhole throat with a minimum radius r(0) = a. 
The range of the radial coordinate is the entire real line, 
— oo < p < oo, and the throat connects two identical 
asymptotically flat space-times. The manifold is every- 
where smooth with scalar curvature given by 



R = - 



2a 2 



(p 2 + ai) 



2\2 ' 



(1.2) 



We consider the problem of computing the self-force on 
a scalar charge at rest in the metric described by ( |1.1[ ), al- 
lowing for the arbitrary values of the coupling constant in 



the wave equation. A similar calculation was considered 
by Khusnutdinov and Bakhmatov [2] who computed the 
electrostatic self-force on a charged particle at rest this 
wormhole space-time. The electrostatic wave equation is 
equivalent to that of a minimally coupled scalar charge 
at rest and hence the self-force on a static scalar charge 
is equal to the electrostatic sef-force, up to an overall 
sign. Bezerra and Khusnutdinov pQ numerically eval- 
uated the self-force on a static scalar for arbitrary cou- 
pling, though we disagree with the overall sign of their re- 
sults. The electrostatic case was reconsidered by Linet [7] 
who derived the Green's function in closed form by trans- 
forming to isotropic coordinates and expanding about the 
Euclidean distance in these coordinates, a method first 
adopted by Copson [5] in deriving the electrostatic po- 
tential in the Schwarzschild space-time. 

The self-force is obtained by taking the gradient of 
the retarded field which is singular at the particle's lo- 
cation and requires regularization. In order to regularize 
the self-force, we compute the Detweiler- Whiting singu- 
lar field, which upon subtraction yields a quantity that 
is regular at the scalar charge's location. We obtain an 
analytic expression for the self-force for arbitrary values 
of the coupling constant, which reveals some expected 
features such as infinite poles which also occurs in the 
expression for the self-force 1] for the wormhole with 
throat profile r(p) = \p\ + a, but we also find some unex- 
pected features such as an infinite set of values for which 
the self-force vanishes. 



II. GREEN'S FUNCTION 

The Green's function satisfies the following inhomoge- 
neous wave equation, 

(□ - £ R)G(x, x') = -g- 1 ' 2 ^ - x'), (2.1) 
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where □ is the d'Alembertian wave operator on the met- 
ric (1.1 1, and £ is the coupling to the scalar curvature 
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which is given by Eq.(1.2). For the charge at rest ac- 



cording to observers moving on integral curves of the 
Killing vector d/dt, the wave equation reduces to the 
three-dimensional Helmholtz equation, 



if 



f 9 ( , 2 2^ 9 \ 1 8 ( . . d \ 1 



(2.2) 



The Green's function may be given as a mode-sum over 
separable solutions to the homogeneous equation, 

G< 3 )(x,x') = ±J2( 2l + V P l( cos 7)9l(p,P% (2.3) 



where Pi(x) is the Legendre polynomial, cos 7 = 
cos 6 cos 9' + sin 8 sin 9' cos A(f) and gi(p,p') satisfies the 
inhomogeneous radial equation, 



dp 



1(1 + 1) 



2£a : 



(p 2 + a 2 ) 



= -5{ P -p'). 

(2.4) 



With the transformation 

V = p/a, 

the radial equation takes a more simple form, 



My-y'). (2.6) 



The general solution may be written as a normalized 
product of two linearly independent solutions of the ho- 
mogeneous equation 



gi(y,y') 



1 M 1] (y<)M 2 \y>) 



N 



(2.7) 



where y< = min{y,y'}, y > = max{ y,y'} and the nor- 
malization constant N is determined by the Wronskian 
of the two solutions. The boundary conditions on the 
Green's function at p — > ±00 require 



^ (1) (y)^0, as 



3 y 



*| 2) (y) -> 0, asy^oo. 



(2.8) 



Writing z — iy, it is clear that the solutions of the ho- 
mogeneous equation (2.6) are associated Legendre func- 



tions of pure imaginary order, 

P^(±iy), Q^(±ty), where p = y/2£. (2.9) 



The multi-valuedness of the associated Legendre func- 
tions gives rise to a discontinuity at y = 0. For y > 0, 
we choose W (y) — Qf(iy) which vanishes as y — > 00 
as required. However, for the solution to be continuous 
across y = 0, we take the branch obtained from the prin- 
cipal branch by encircling the branch point z = 1 (but 
not the point z — —1) once. If we denote this branch by 
Qi^z), then it can be shown [5] 



Then, the function 



(2.10) 



y>0, 



*| 2) (2/) = <! 



r(Z - M + 1) 



y<o, 



is continuous for — 00 < y < 00 and satisfies the ap- 
propriate boundary condition and hence is the correct 
choice. 

The symmetry of the space-time implies that we take 
our inner solution to be (y) = ^i 2 \—y), or given 



(2-5) explicitly 



M 1) (y) = { 



y>o, 

y<0. 



This solution clearly satisfies the vanishing boundary 
condition at y — > —00 (p —¥ —00) and is linearly indepen- 
dent to iffj (y) in the two regions of the space-time. We 
can re- write the solution in a more symmetric form using 
standard results for the Legendre functions [TU] , yielding 



^\y) 



(-!)<+ V™ Q?(iy)+™ 



T{l + p + l) 
T{l-p + l) 



y>o, 

y<0. 



The normalization constant is given by 

N = -(y 2 + l)W{¥ l 1) (y)^\ 2 \y)} 

= 7r( _ iy+le2 ^r|±M±il. 

T(l - p + 1) 



(2-11) 
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So the radial solution may be written as 

gi(y,y') 

a 7r a 



where 



a ir 



-Qr"(«/)Qf(%') 



if y,y'>o, 



a 7T 



: Qr^(iy)Qr(ij/o 



if v, y' < o, 
if y > 0, y' < 

or y' > 0, y < 0. 



Let us assume, without loss of generality, that y, y' > 
(the closed-form expression won't depend on this condi- 
tion), then the Green's function is 

oo 

G( 3 >(x,x')= *£)(2Z + 1)^(0087) 



Z=0 



e-^'p,-"^)^ (%>) + - e ^gr M (^)or(^') 

7T7 



(2.12) 



In a recent paper , we have derived the following sum- 
mation formula for associated Legendre functions of ar- 
bitrary complex order, 

oo — /icosh _1 (x) 

+ l)P(cos7)e-^Pf "(r^Qf {v>) = ^— ^ 

z=o 

(2.13) 

where 



R = rf + rj' 2 — 27777' cos 7 — sin 2 7, 



X 



rjr] — cos 7 

(77 2 - 1)1/2(7/2 _ 2)1/2 



(2.14) 



and where the radial variable here is real and runs over 
the range 77 > 1. We can analytically continue this result 
by making the transformation 77 = iy and taking the 
appropriate branch. We obtain 

00 

X> + i)p ; (cos7) e -^ ,; pr Al (^<)Qr(^>) 
1=0 



' g/^cosh 1 (x) 

1R 1 / 2 

e -Mcosh _1 (x) 
iW 2 

g-/«H e /iCOsh _1 (x) 

iR}T 2 



y,y'>o, 
y,y' < 0, 

y > 0, y 1 < or y < 0, y 1 > 0, 



R = y 2 + y' 2 - 2yy' cos 7 + sin 2 7, 
yy' + cos 7 

X = ( y 2 + 1 )l/2(^2 + 1 )l/2- 



(2.15) 



These results may be checked numerically by multiplying 
both sides by P//(cos7) and integrating with respect to 
7- 

For y,y' > 0, we can use the well-known relations be- 
tween Legendre functions to rewrite the Green's function 
| |2.12| ) as 



G( 3 )(x,x') 



87r|a| sin(/i7r) 



£(2Z + l)P,(cos 7 ) 



1=0 



- 2 ^P^iiy<)Qt(iy>) + e 2 ^pnw<)Qr{w>) 

(2.16) 



We can employ the summation formula (2.151 to obtain 



the following closed form representation of the Green's 
function for a static scalar charge in our wormhole space- 
time 

sin(^cos _1 (-x)) 



G( 3 )(x,x') = 



1 



(2.17) 
(2.18) 



47r|a| sin(/j7r) R 1 / 2 

where we have used the fact that 

7r + i cosh - (x) = cos _1 (— x) f° r \x\ < !• 

A similar analysis for the case where y,y' < and where 
the two points are in different regions of the space-time 
yield the same closed form expression as Eq.(2.17l. Fi- 
nally, restoring the variable p = y/a gives 

G( 3 >(x,x') 

1 sin(/zcos _1 (— x)) 

An sm(pn) (p 2 + p' 2 - 2pp' cos 7 + a 2 sin 2 7) 1 / 2 ' 

(2.19) 



where 



X 



pp' + a 2 cos 7 



(2.20) 



(p 2 +a 2 y/ 2 {p 12 + a 2 ) 1 /2' 

We note that this Green's function has the correct 
Hadamard singularity structure |12j in three dimensions. 
We also note that the Green's function possesses infinite 
poles for certain values of the coupling constant given by 

n 2 

£=— , neZ\{0}. (2.21) 

This is analogous to the infinite poles that arise in the 
Green's function for throat profile r(p) — \p\+a, as shown 
in Ref. pQ. 



For minimal coupling, Eq.( 2.19 ) is understood to mean 

G (3) = 1 COS-^-x) 

M=0 4tt 2 (p 2 + p' 2 - 2pp' cos 7 + a 2 sin 2 7)1/2' 

(2.22) 

which agrees with the closed form representation given 
in Ref. [5] if we take 7 = 0. 
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III. SELF FORCE AND SELF ENERGY 

There are a number of approaches one can take to 
compute the self-force (see [13] for example), but the 
most direct way for a static charge is by subtracting the 
Detweiler- Whiting [Hj parametrix which yields a finite 
quantity upon taking coincidence limits. It is constructed 
in such a way as to leave, upon subtraction from the 
Green's function, only the radiative part of the field en- 
tirely responsible for the self-force. This approach was 
adopted in Ref. [15] to compute the self-force on a static 
scalar charge in Kerr space-time. There it was shown that 
for a static scalar charge with world- line x' a = z a (r) and 
four- velocity u a = u 5 a t> in a general stationary space- 
time, the self-force is 



/- lf = 47rg 2 lim [v Q ( — (x, x') - G DW (x, < 



(3.1) 

where G^ 3 -*(x, x') is the zero-frequency mode of the four- 
dimensional retarded Green's function (modulo a factor 
of 2ir) and Gnw(x,x') is the Detweiler- Whiting Green's 
function given by [T3] 

G DW (X X') = 1 / A^S.Srct) + AV^ar.gadv) 

47r V 2 r, fit 2 r ar iv 



V{x, X '{T))dT 



where A is the Van-Vleck Morrette determinant, r rot is 
the retarded distance between the field point x and the 
retarded point x 1 = x Te t on the charge's world-line and 
?"adv is the advanced distance between x and the advanced 
point x 1 = :r a dv on the world-line. 

In the case of minimal coupling, it has recently been 
shown [16] that the Detweiler- Whiting Green's function 
for a static charge in a static space-time is equivalent 
to the direct part of the Hadamard Green's function on 
the spatial part of the metric up to the order required 
for regularization, and equal up to all orders for ultra- 
static space-times such as the wormhole space-time under 
consideration. In that case, the appropriate parametrix 
that must be subtracted to give the correct self-force is 
particularly simple, 



G D w(x,x') 



1 Ai/^x.x') 



(3.3) 



where A and the world-function a are calculated on the 
three-dimensional metric 



ds 2 



(3) 



= dp 2 



\ P )dn 2 



(3.4) 



However, it is clear that this statement cannot be true for 
arbitrar y co upling since the biscalar V(x, x') appearing 
in Eq. (3.2 1 depends on the coupling constant £ while 



the direct part of the three-dimensional Hadamard form 



(3.3) does not, indeed it is purely geometrical. 



In obtaining an expression for Gdw(s, x'), we will use 
the standard expansions 

A 1 ^ 2 (x,x') = 1 + ±R ab a> a a> b + 0(a 3 / 2 ) 

V(x,x') = ^-i)R + 0(a^ 2 ). (3.5) 

We can use the fact that the space-time is spherically 
symmetric to set 7 = 0, in which case the Van-Vleck 
Morrette expansion reduces to 

Al/ W) = l-^^Ap 2 +<W). (3-6) 



From the expansion for V(x,x'), the integral in Eq.(3.2) 



V(x,x'(r))dT= |(e-i)i?Ar + 0(Ax 2 ), (3.7) 



where At = T a( jv — Tret- Coordinate expansions for the 
quantities r rct , r a d v and At were computed for a static 
charge in an arbitrary stationary space-time in Ref. [15] . 
which in the wormhole space-time under consideration, 
simplify greatly to give 



/-,-, 1 : / a ,iv = I At = \Ap\. 
Hence the Detweiler- Whiting parametrix is 



(3.8) 



(3.2) G DW (p,p') = 



1 p 2 a 2 



4n\Ap\ (p' 2 + a 2 y 



:\Ap\+0(Ap 2 ), (3.9) 



where recall that p — y/2£ and we have used the expres- 
sion for the scalar curvature given by Eq.(1.2). 



We also require an expansion for the closed form 
Green's function of Eq.(2.19). Taking the scalar parti- 



cle's location to be p' and the field point to be p, then 
for p near p' the Green's function may be expanded as 



G( 3 W) 



1 



1 p COs(/17r) 



AirAp 4-7T sin(/i7r) (p' 2 + a 2 ) 
1 \a\p /2p' cos(/x7r) + |a|/isin(/i7r) 



8ir sin(/i7r) V 
0(Ap 2 ), 



(p> 2 + a 2 ) 2 



Ap 
(3.10) 



where Ap = p' — p which we assume, without loss of 
generality, to be positive. 

The self-energy for a charge q is proportional to the 
coincidence limit of the radiative field, and is given by 



U(p') 



2nq lim 

q 2 p cos(pir) 



G^(p,p')-G B w(p,p') 



2 sin(^7r) (p' 2 + a 2 )' 



(3.11) 



For minimal coupling, £ = (p, = 0), we arrive at the 
result of Khusnutdinov and Bahkmatov [5] (which was 
rederived by Linet [7]), 



U(p') 



q 2 \a\ 



2TT(p' 2 + a 2 )' 



(3.12) 
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FIG. 1. Plot of the radial component of the self-force on a static scalar charge in a wormhole space-time. We have set q — 1 
and o = l. The thickness of the curves increases as the coupling increases and the continuous lines are for values in the range 
€ < 1/2- 



while for conformal coupling £ = 1/8 (// = 1/2), the self- 
energy vanishes, which agrees with the results of Bezerra 
and Khusnutdinov 1]. The self-energy, like the Green's 
function, is also divergent for the particular values of the 
coupling constant given by Eq. (2.21). 



Turning now to the calculation of the self-force, which 
we have defined for a static scalar charge in a stationary 
space-time in Eq.(3.1). For the ultra-static wormhole 



space-time under consideration, this reduces to 
f eU = W lim V p (g<® (p, pi) - G dw (p, P r 



(3.13) 



with all other components vanishing due to the spherical 
symmetry. Substituting in the expression for the singular 
field (3.9| and the Green's function expansion (3.10), we 



obtain for the self-force 



/; clt = g>|MCot(/i7r) 



P 



(p l2 + a 2 ) 2 ' 



(3.14) 



Again, we note that the self-force has an infinite number 
of poles whenever the coupling constant is £ = n 2 /2, 
where n € Z/{0}. A similar analysis [1] for throat profile 
r (p) = \p\ + a a l so exhibits this behaviour for the self- 
force. 

For minimal coupling, we obtain 



j^self 

Jp — 



q 2 p'\a\ 



n(p'' 



(3.15) 



which is in agreement with the electrostatic self-force 
derived in Ref. [2], modulo the sign of the force since 
the field for a minimally coupled static scalar in an 
ultra-static space-time is minus the electrostatic field. 
Hence the self-force on a static scalar in the wormhole 
space-time is always repulsive with respect to the throat 
whereas in the electrostatic case, it is always attractive 
with respect to the throat. 

For the conformal coupling in three dimensions , £ = 
1/8, we obtain zero self-force which is in agreement with 



the result of Ref. pQ. In fact, just as there are an infinite 
set of values for which the self-force diverges, there are 
also an infinite set for which it vanishes, 



/f lf =0, for/i = 



2n + 1 



(2n + l) 2 



(3.16) 



That the self-force vanishes for £ = 1/8 for a massless 
field is expected since the spatial section of the metric is 
conformally flat. However, the existence of an infinite set 
of values of the coupling constant for which the self-force 
vanishes is a surprising result, at least to this author. 

In Fig. [TJ we plot the self-force for various values of 
the coupling constant. The first thing to note is that 
the overall sign of the self-force differs to the numerical 
computation of Ref. [1 . So for the field with conformal 
coupling in the range < £ < 1/8, the self- force is every- 
where repulsive with respect to the wormhole throat, not 
attractive as previously claimed. The self-force vanishes 
for conformal coupling £ = 1/8 and then becomes attrac- 
tive for values in the range 1/8 < £ < 1/2. The force be- 
comes increasingly attractive with respect to the throat 
as £ increases towards £ = 1/2 where it becomes diver- 
gent. This cycle then continues, the force is repulsive for 
1/2 < £ < 9/8, vanishing at £ = 9/8 and attractive with 
respect to the throat for 9/8 < £ < 2 and diverges for 
£ = 2 etc. 

These poles in the expression for the self-force, which 
also appear in the closed form representation of the 
Green's function, can be understood in the context of 
quantum mechanical scattering theory as the energy of 
bound states. In this classical context however, the ex- 
istence of such poles suggests that for certain coupling 
strengths of the scalar field to the curvature, the charge 
cannot remain static. 

For a particular value of the coupling constant, the 
direction of the self-force in either part of the space-time 
is independent of the radius at which the scalar charge is 
being held. For example, for 1/8 < £ < 1/2 the force is 
always attractive with respect to the throat regardless of 
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where in the wormhole space-time the charge is placed, 
and hence for the coupling strength in this range we may 
hypothesize that the scalar charges will accumulate in the 
vicinity of the throat. On the other hand, for £ < 1/8, 
the self-force is everywhere repulsive with respect to the 
throat. The magnitude of the self-force for a static scalar 
is maximized when the charge is held at p = ±|a|/\/3, 
and the force at this radius is 

.self _ i 3^3gVcot(^7r) 

/max-± " 16fl2 • {^V 



IV. CONCLUSIONS 

We have obtained an analytic expression for the 
self-force on an arbitrarily coupled static scalar charge 
in a wormhole space-time with throat profile r(p) = 



p 2 + a 2 . Analytic expressions have previously been ob- 
tained only for minimally and conformally coupled scalar 
fields, while the case of general coupling had been com- 
puted numerically. We find that there are infinite poles 
in the expression for the self-force corresponding to the 
values of the coupling constant where £ = n 2 /2, for 
n e Z \ {0}. We also find that there are an infinite set 
of values of the coupling constant for which the self-force 
vanishes. 
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